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Abstract 

In this note we investigate bound states, where scalar and vector bosons are trapped by BPS vortices in 
the Abelian Higgs model with a critical ratio of the couplings. A class of internal modes of fluctuation 
around cylindrically symmetric BPS vortices is characterized mathematically, analysing the spectrum of the 
second-order fluctuation operator when the Higgs and vector boson masses are equal. A few of these bound 
states with low values of quantized magnetic flux are described fully, and their main properties are discussed. 
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1. Introduction 

Very soon after the discovery of Abrikosov quantized flux lines in the Ginzburg-Landau theory of Type 
H superconductors [T], the existence and nature of fermionic bound states on these vortex filaments were 
discussed by de Gennes et al. in Reference [ 2 ] . Quantized magnetic flux lines were rediscovered by Nielsen 
and Olesen in the Abelian Higgs model, see [ 3 ], a finding that enhanced the interest of these topological 
defects by promoting them to the relativistic and quantum world. By adjusting the couplings in the Abelian 
Higgs model to drive the system to the critical point between Type H and Type I superconductors, Bo- 
gomolny showed, see [ 3 ], that quantized vortex lines still exist but move without interaction with respect 
to each other. Bosonic vortex bound states were investigated by Goodman and Hindmarsh, see [S], in the 
context of the Abelian Higgs model for any value of the parameter governing the transition between Type 
I and Type H superconductivity. In this short note we shall focus on finding BPS vortex bound states and 
we shall describe these internal boson-vortex modes by a mixture of analytical and numerical methods, at 
least at the same level of numerical precision as the BPS vortex solutions themselves. 


2. BPS vortex fluctuations 


The Abelian Higgs model describes the minimal coupling between a C/(l)-gauge field and a charged 
scalar field in a phase where the gauge symmetry is broken spontaneously . In fact, it is a relativistic 
version of the Ginzburg-Landau theory of superconductivity. At the transition point between Type I and H 
superconductivity, where the masses of the Higgs and vector fields are equal, the AHM action reads 







( 1 ) 


Here, non-dimensional coordinates, couplings and fields are used, while (p{x) = ^i(a;) -I-1^2(2;) is a complex 
scalar field and A^{x) = {Aq{x), Ai{x), A2{x), A^{x)) is the vector potential. The covariant derivative is 
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defined in the conventional form, Dfj_(j){x) = {d^ — iA^{x))4>{x), whereas the electromagnetic field tensor is 
also standard: F^i,{x) = d^Aij{x) — dvA^{x). We choose the metric tensor in Minkowski space as = 
diag(l, —1, —1, —1), = 0,1,2, 3, and use the Einstein repeated index convention. In the simultaneous 

temporal and axial gauges = ^3 = 0, the Bogomolny arrangement of the energy per unit length for static 
and xs-independent field configurations E[i 5 i), ^], see [1], shows that solutions of the first order PDE system 

Di(j)±iD2^ = 0 , Fi2±i(</)>-l)=0, (2) 

with appropriate asymptotic behaviour at infinity in the Xi : a: 2 -plane, are absolute minima of A^. It 
was proved in [5] that there exist solutions of the PDE’s with finite string tensions that are proportional 
to the magnetic flux along the xa-axis of n G Z quanta: E[(/),A] = ^ | (PxFi 2 \ = TT\n\. These topo¬ 
logical objects are denoted BPS vortices because they correspond to the Abrikosov-Nielsen-Olesen vortex 
filaments arising in Type II superconductors, see m-m, when the scalar and vector penetration lengths in 
the Ginzburg-Landau free energy are equal and the system lives exactly at the transition point to Type I 
materials. 

Denoting the BPS vortex fields as: 


(j)v = i^ix;!!,) = tpi{x-,n) + iip 2 {x]n) , Ay = {Vi{x;n),V 2 {x-,n)) with x = {xi,X 2 ) , 


and assembling the vector and scalar vortex fluctuations (ai(x), a 2 (x)), (p(x) = ‘pi(x) +i(p 2 (x) in a column 
vector ^(x) with transpose 


^*(x) = ( ai(x) 02 (x) <pi(x) (p 2 (x) y , 


one checks that the linearized dynamics is governed by the action of the second-order vortex fluctuation 
operator 


( -A + IV'P 
0 

—2Di'ip2 

V 2Di^i 


0 -2Ditjj2 \ 

-A -b IV'p -2021^2 2D2V'i 

-2D2i^2 -A+i(3|V^|2-l) + I4I4 -2Vkdk-dkVk 

202^1 2Vkdk+dkVk -A+i(3|^|2-l) + EfcEfe / 


(3) 


on ^{x). Resolution of the spectral problem 'H'^^x{x) = uj‘l^\{x), where A is a label in either the discrete 
or the continuous spectrum of , permits the decomposition of ^(af) as a linear combination of the ^\{x) 
eigenfunctions. 

In the search for bound state (normalizable) eigenfunctions other than zero modes, i.e. 0 < < 1 

assuming that = 1 is the scattering threshold, we shall profit from a hidden SUSY structure of 
Linear deformation of the PDE system ([^, together with the background gauge 

'9lK(^)]l +92[C(S“)]2 - Iplix) [^{x)]4, + 1p2{x) [^{x)]3 = SfeOfc (f) - -01 (x) </52 (:?) + V'2 (A) <^>1 (£) = 0 , (4) 


is encoded in the following first-order PDE operator 
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acting on the space of BPS vortex fluctuations. This operator V allows us to embedd in a SUSY 
Quantum Mechanical system because = 'D^'D and we find the following SUSY algebra generated by the 

supercharge Q = 


0 0 
V 0 


2 






% being the SUSY Hamiltonian, while the SUSY partner to is: 

A + |V^|2 0 0 0 \ 

0 -A+|V^|2 0 0 

0 0 -A + i(|^|2 + l) + UfeUfc -2Vkdk-dkVk ■ 

0 0 2VA + dkVk -A+i(|V'|2 + l) + UfcUfe / 

Except for the eigenfunctions in the kernel of V, which are zero modes of T-L^ ^ the two operators are 
isospectral. This supersymmetric structure led to the proof of the Weinberg index theorem on the plane 0 : 

f '1 

ind V = dim Ker T) = lim Tr ^2 < -—— — ——r-—-r > = 2n , 

M^oo ^ I X>tx» + M2 x>x>t + M2 J 

stating that has 2n zero modes in its spectrum: = 0, I = 1,2, - ■ ■ ,2n. Moreover, ^ also 

guarantees that the 'H'^-spectrum is non-negative, such that the zero modes of V are all the ground states 
of H because KerH“ = 0. 

We shall focus on BPS cylindrically symmetric vortex filaments shaped according to the Nielsen-Olesen 
ansatz: 

</'(^) =/n(r)e“® ; rAg{r,9) = nPnir) . (8) 

We stress that: (1) Cylindrical coordinates are chosen in the K^-space and the vector potential components 
are adapted to them. (2) Besides the temporal and axial gauges, the radial gauge = 0 is assumed such 
that the vector field is purely vorticial. (3) The complex scalar field is expressed in polar form. 

The first-order PDE system ([^ becomes the following ODE system: 

= ^/„(r)[l-/3n(r)] , ^(r) = ^[l-/2(r)]. (9) 

The solutions for the radial profiles /n(r) and /?«(?'), in the xi : a: 2 -plane and infinitely repeated along the 
cca-axis, determine the cylindrically symmetric BPS vortex solutions. The finiteness of the energy per unit 
length demands that fn{r) —t 1 and /3„(r) —)■ 1 as r —)■ oo. 

Some analytical progress in the investigation of the zero-mode fluctuations on BPS cylindrically sym¬ 
metric vortices was achieved in [7]. Further comprehension of their structure was obtained in References 
[5] and [5]. Here, the motivation leading several researchers to describe in detail the vortex zero modes 
came from the study of vortex scattering at low energies within the approach of geodesic dynamics in their 
moduli space, see e.g. m- In this note we shall focus on finding and describing excited fluctuation modes in 
the discrete "H^-spectrum, i.e., internal modes of fluctuation, where the BPS vortex captures scalar and/or 
vector mesons, an issue not discussed in the literature on the Abelian Higgs model. 

3. Spectrum of cylindrically symmetric BPS vortex fluctuations 

In the search for positive bound states Ca (^) discrete spectrum of the operator , the use 

of supersymmetry is convenient. If > 0, the SUSY structure M implies that 

Moreover, the eigenfunctions of are related through the supercharges: {x) = ^ (x). In addition, 

it should be recalled that the background gauge condition must be satisfied in order to eliminate spurious 
gauge fluctuations. The strategy is thus to solve the spectral problem for 'H~ first and apply the 
operator to the eigenfunctions, finally obtaining the eigenfunctions of . This indirect path is more 
appropriate because the spectral problem of 'H~ is more tractable owing to its block-diagonal form. Two 
classes of eigenmodes of the operator 'H~ can be distinguished: 

• Class A -eigenmodes: The two 1x1 block-diagonal sub-matrix differential operators in 'H~ prompt a 
complete decoupling of the vector field from the scalar fluctuations in the ^“-spectral problem. There thus 
exist eigenfunctions of the form = ( ai{x) 0 0 0 )* and [Ca~(^)]* = ( ^ a 2 {x) 0 0 )*. 

The H^-wave functions of the form do not satisfy the background gauge. Therefore, we shall 


(- 

= 

V 
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study only the physically meaningful possibilities (a;)] among this class of eigenmodes. The non-null 


component of (a;) complies with the PDE (—-1- \4’\'^) ai(al) = oj\ ai{x), or, in polar coordinates, 


dr'i 


1 dai 
r dr 


1 d'^ai 

j,2 Q()2 


+ ifnir) - ujl]ai =0 . 


( 10 ) 


The separation ansatz ai(a;) = Vnkij") cos{k0) [pleads to the ID Sturm-Liouville problem 


_ <Pvnk{r) _ 1 dvnkjr) 
dr^ r dr 

for the radial form factor u„fc(r). Univaluedness of the fluctuations demand that the wave number k must be 
a natural number: A: = 0,1, 2 .... The ODE (11) is no more than a radial Schrodinger differential equation 
with a potential well V^{r;n,k) = f^ij) + which includes a centrifugal barrier when k ^ 0, bounded 
below and running to 1 at infinity: limr_).oo = 1. Consequently, a continuous spectrum arises in the 

ujf G [l,oo) range, i.e., for energies above the scattering threshold = 1. Below this threshold, in the 




fn(r)-^^x + Vr,k{r) = 0 


( 11 ) 


uj\ G (0,1) range, boson-vortex bound states may exist if the spectral problem (11) admits eigenvalues. The 


procedure to find both the eigenvalues a;| and the eigenfunctions u„fc;j(r) will be implemented in the next 
Section for low values of n. The need to identify the eigenvalues will leads us to convert the ODE 0’ 
where is a priori unknown, in a system of equations of finite differences by some discretization method 
of the half-line to a lattice with a finite but large number of points. Diagonalization of the matrix of the 
linear system in turn provides the eigenvalues and eigenfunctions, which are very good approximations to 
the eigenfunctions and eigenvalues of 'H~ provided that the number of points of the discretization is large 
enough. 


Class B "H -eigenmodes: The 2x2 block-diagonal sub-matrix in Ti. 

(p 2 (x) y, leading to the spectral PDE system: 


acts only on scalar field fluctuations 


of the form 


B- 


(f)]‘ = ( 0 

[-V^T 


0 


1)+ VkVk](pi{x)-2Vkdk(p2{x) = ujlipi{x) 
l) + VkVk](p2{x) + 2Vkdk(pi{x) = u}lip2{x) . 


( 12 ) 

(13) 


Eor cylindrically symmetric BPS vortices, the ansatz (pi(x) = — Unk{x) sin[(fc-|-1)0], '^ 2 {x) = u„fc(r) cos[(fc- 


1)0] that converts the PDE system (121 and (13) into the spectral ODE is clear: 


Tlunkir) = - 


d^Unk{r) 1 dunk , fl 


L2(l+/n(^)) + 


in/3n{r) - (1 -k fc))2 


ujf Unk{r) = 0 


(14) 


dr^ r dr 

where the radial form factor Unk{x) is the unknown. This is a radial Schrodinger differential equation with 
an effective potential well: Vgg(r) = i[l -k /^(r)] -k :^[n/3nir) — (1 + k)]'^. Erom the functional behavior of 
Vgg, we may conclude that a continuous spectrum arises settled on the threshold value oj^ = 1. The same 
reasoning indicates that there could exist bound states in this class with eigenvalues in the uj G (0,1) range. 
However, a theoretical ar gum ent can be used to discard this possibility. The linear differential operator TZ 
associated with the ODE (14) can be factorized as 7^ — 1 = where £ refers to the first-order differential 


operators C = —^-k^[l-kA: — n/3„(r)]. This means that the 7^—1 operator has a non-negative spectrum and, 
consequently, there are no bound states in the discrete spectrum of 'H~ within this class B of fluctuations. 

Translation via use of the supercharges of all the 7A~-spectral information described previously reveals 
the structure of the TA'^-spectrum: 


• Class A eigenmodes: Assuming knowledge of Lof and Vnk{r) from the solution of (11) the eigenfunctions 
of paired through supersymmetry with these class A 7A~-eigenmodes take the form: 

/ sin0cos(fc0) ^""^^^’~^ - ^Unfc(r)cos0sin(fc0) \ 




A-l- 


(a:, n, k) = 


— cos 0 cos(fc0) — yu„fc(r) sin0sin(A:0) 

fn{r)v„k{r) cos(n9) cos(fc0) 
fn{r)vnk{r) sin(n0) cos{k9) 


fc = 0,l,2,... . 


( 15 ) 


^ ai{x) = Vnki'f') sin(A:0), k = 1,2..., also leads to 0- We shall pursue the cosine alternative, for the sake of brevity. 
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It is immediate to check that satisfies the background gauge Q, meaning that the fluctuations (151 

correspond to admissible eigenfunctions of the second-order BPS vortex fluctuation operator 'H'*' From 
(111 we may conclude that a continuous spectrum emerges at the threshold value = 1, while the discrete 
spectrum is confined to the open interval ujj G (0,1). The presence of bound states in the ^'''-spectrum, 
however, requires the existence of eigenfunctions of ( |ll[ ) satisfying the boundary conditions ^^^(0) = 0 and 
limr_>.oo Vnk{r) = 0. This point will be addressed in tEe next section. 


• Class B H+-eigenmodes: The corresponding SUSY partner 'H’*'-eigenfunctions associated with the 

class B ^“-eigenfunctions are given by: 


= r” 


/ hnkir) sm[{n - k - 1)9] 
hnk{r) cos[(n- k- 1)0] 
cos{k0) 


V 


U{r) 
■7^sin(fc0) 


\ 


The new radial form factor hnkif) may be defined from the relation M„fe(r) = hnk{f) in such a way 

that the ODE (14) turns into the equation 

r Kk(r) + [-l-2k + 2n/3„(r)]/i;fc(r) -f - /^(r)]/i„fc(r) = 0 , 


for the radial form factor hnkir). Regularity at the origin, however, of the positive eigenfluctuations 
requires that 0 < k < n — 1. In this case, there exists a continuous spectrum emerging from the threshold 
value = 1 and no positive eigenvalue bound states arise. Notice, however, that the form of these 
eigenmodes follow the ansatz given in [5] for the zero modes, such that we can regard the 2n zero 

modes as the only bound class B H'^-eigenmodes[^ 

Orthogonality between eigenfunctions belonging to different classes is guaranteed by the conservation of 
the scalar product in the SUSY partnership: (x), [x)) = (x), i^)) — because class A 

and B eigenfunctions of 'H~ are clearly orthogonal. Orthogonality between eigenfunctions belonging to the 
same class with different angular dependence is established by Fourier analysis. 


4. Positive eigenvalue bound states of the small vortex fluctuation operator 


We now attempt to elucidate the existence of excited fluctuations of class A belonging to the discrete 
spectrum of T-L^ with positive eigenvalues lower than 1. The search for and the analysis of these fluctuations 
reduce to the numerical computation of the radial form factor Vnk{r) in the ODE Our strategy to 

achieve this is to employ a second-order finite-difference scheme that simulates the differential equation 0 
by the recurrence relations 


^nk\j ^^nk\j 


(i-l) 

^nk-,j 




-^nk\j 


- V, 


nk;j 


(Ax)2 


2f(Aa:)^ 


+ 


fnii^x) 


'(Aa;)2J 


^nk;j 


_ 2 i'i) 

^nk;j ^nk;j 


(16) 


where we have confined the problem to the interval [0,rmax] for a large enough r^ax- We denote = 
Vnk-jii^x), with Ax = and choose a mesh of N points with i = 0,1,2, ■■■, N. The eigenfunctions and 
the eigenvalues depend on the values of the angular momentum k and the vorticity n. The index j is used 
to enumerate the discrete eigenfunctions. The contour conditions are: 


( 1 ) - 


4 V. 


( 2 ) 

nk\j 


— V, 


(1) 

nk\j 


(Ax)2 


fli^x) + 


(Aa 


( 1 ) 2 ( 1 ) 
^nk-j = ^nk-j 


and (2) 


^nk;j 


= 0 


^The sine alternative leads to degenerate ei genf unctions with {x,n, k) obtained by simply replacing cos(k9) by sin(A:0) 
and sin(/c0) by — cos{k6) in the wave function |T^. Knowledge of the radial form factor Vnki't') authomatically gives both the 
cos and sin bound states. 

^Other n zero modes are easily generated by rotating 7r/2 separately in the scalar and vector field fluctuations, which is a 
symmetry of the spectral problem. The same argument can be applied to the eigenfunctions 
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A good estimation of the discrete eigenvalues is obtained through diagonalization of the N x N matrix 
in the left member of the linear system (fb). We show the eigenvalues of "H”*’ for low values of n and k 
obtained in a Mathematica environment in Table 1 by applying this procedure with the choice of = 400. 
In Figure 1 we illustrate the behavior of the potential wells of the radial Schrodinger equation 0 for n = 4 
and the overlapped dashed lines that determine the discrete eigenvalues of the 'H~*'-spectrum in this case. 
The radial form factor Vnk;j{T^) of the eigenfunctions associated with these eigenvalues are also shown. 

In general, we observe that the number of bound states increases with the magnetic flux n. In particular, 
we conclude the existence of I bound state for n = 1-vortices; 2 bound states for the n = 2 and n = 3 
vortices; 3 bound states for n = 4-vortices, and 4 bound states in the case of n = 5-vortices. In Table 1 we 
include a graphical representation of the discrete T^’^-spectrum for several values of the vorticity n. 


Eigenvalues of the discrete spectrum of 

n 

O 

II 

k = 1 

k = 2 

1 

^ 0.777446 

- 

- 

2 

(wyi)2 = 0.538573 

(a;A.^)2 = 0.972563 

- 

3 

(a;A,^)2 = 0.402692 

(a;A.^)2 = 0.830078 

- 

4 

^ 0.319276 
l^to-,2? = 0.988212 

(a;A.^)2 = 0.701708 

- 

5 

(cjA ^)2 = 0.263671 
(cuA I 2 )" = 0.939461 

(a;A,^)2 = 0.601223 

(w4;i)^ = 0.942438 


Table 1: Numerical estimation of the discrete spectrum eigenvalues associated with the class A eigenfunctions n, k) with 

angular momentum k, together with a graphical representation of the second-order small n-vortex fluctuation operator 
spectrum. The eigenvalues of the form cOnO shown in this Table agree with those given in Table 3 of Reference for critical 
quotient of the scalar and vector field masses up to a factor of 2 due to a different convention. 




Figure 1: Plots of the effective radial potential wells V"g^(r, n = 4, /c) arising in Jill for k = 0,1,2 (solid lines) with the discrete 
eigenvalues (dashed lines) overlapped (left) and the corresponding radial eigentunctions (right) of the small n = 4 

vortex fluctuation operator . 


In sum, there exist bound states eigenfunctions of . Here, we have described the 

stationary wave functions where an awkward combination of scalar and vector boson fluctuations are trapped 
by a cylindrically symmetric BPS vortex. These configurations oscillate in time with frequencies determined 
by the discrete eigenvalues and are thus internal modes of fluctuation of the BPS n-vortex. 
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